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We study the effect of interactions on the time reversal invariant topological insulators in four 
and three spatial dimensions. Their topological indices are expressed by the interacting Green's 
functions. Under the local self-energy approximation, we find that interaction could induce non- 
trivial frequency-domain winding numbers and change the topological classes of the system. Our 
results suggest that the topological phases could be destroyed without developing long range or- 
ders. Practical issues on the accurate frequency-momentum integration combined with DMFT and 
diagrammatic calculations of the interacting Green's functions are also addressed. 
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I. INTRODUCTION 



OQ Rec ently, there has been a raising interest on a class of topological phase of matter: the topological insulators, 

see^S for reviews. Questions like what is the interaction effect on the topological insulators and how to characterize 

I I ' the topological phase beyond single particle basis pose intriguing challenges. There are some previous studies on the 

interacting effect on topological insulators. Wang et al."^ performed exact diagonalization and quantum Monte Carlo 
^ studies of the interacting Haldane model. Hohenadler et al.Hl, Zheng et al.^^, Yu et alM performed quantum Monte 

"t^ Carlo and variational cluster approach studies on the Kane-Mele-Hubbard (KMH) model. One interesting finding in 

1 ; these studies is that the topological phases give their way to a gapped featureless state (spin liquid phase for KMH 
Cd model) at small value of the spin orbital coupling. The possibility of fluctuation effect driven by interaction breaks 
H the topological phase significantly go beyond the mean-field treatments''^. Since in the latter case, long range orders 

I (LRO) are needed to compete with the topological phases. 

However, one unsatisfactory feature in these studies is that there lacks a direct characterization of the topological 
order: the phase boundaries are mostly determined by indirect signatures such as the correlation functions for LRO 
fj or the excitation gaps, rather than a sharp change of a topological indices. This leaves the mechanism and the nature 
of the fluctuation induced topological transitions unclear. 

In this paper, we study the interaction effects on the topological insulators based on their single-particle Green's 
functions^. We start fromthe four dimensional Chern insulator^ since it is considered as the mother state of large class 
-y-K of topological insulators® ^°. The system's topological index- Chern number- is expressed by the following formula: 

o 

^. n= ^^Tv f ^ r t^ Gd,G-'Gd.G-'Gd,G-'Gd,G-'GdrG-' (1) 

I where Sfj^upaT is the five-order anti-symmetric tensor, pl,u...t indices denote the frequency- momentum 

^-H (w, kx, ky, kz, k\) indices. The frequency integration ranges from — oo to oo and the momentum integration is over the 

first Brillouin zone. 

To the authors knowledge, in the quantum Hall effect context the formula was first appea red irPilandi^l. It was latter 

S^ used to study the finite temperature and disorder effect on integer quantum hall effedP^ini. It was also used extensively 

?-H by Volovik in his book^^. Recently, Ref.® ^'^ proposed to use it and its variants as the topological order para meters 

^ of the interacting topological insulators. It is also used to study the boundary effect of topological insulators by^^l^. 

Physically, the formula descries Hall conductance of the material since it appears as the coefficient of the Chern-Simons 

term. On the mathematical side, the outcome of the integration EqIT]is guaranteed to be an integer with physical 

Green's functions for a gapped system since it represents the homotopy index tt^{GL{M, C)) = Z. 

For noninteracting models G~^(w, k) = iuj — H]^. The frequency integration could be done analytically, results 
into the Chern number^'' expressed by momentum space integrations, i.e. the TKNN numberglSI. For interacting 
case, the Green's function becomes G~^{uj, k) = iw — H^ — S(a;, k). The self-energy T,{uj, k) contains the information 
of many-body physics. The fluctuation effect of many body physics are encoded in the frequency dependence of 
the self-energy. In general, the five-fold frequency-momentum integration in EqnlT] becomes involved and one can 
not complete the frequency integration analytically. In mean- field studies of the interaction models ■ and Born- 
approximation study of the disorder effects on topological insulators^fl', the self-energy is approximately frequency 



independent (S(a;,k) w Sk) and it is effect is effectively absorbed into i?k and modified the topological classes 
thereof. However, more ubiquitous and intriguing cases are the self-energy has strong frequency dependences where 
the many body effect does not simply renormalize the single particle Hamiltonians. The main goal of the present 
paper is to study how nontrivial frequency dependency of the self-energy affects the topological classes. In section 
ITT] we introduce the local self-energy approximation and study its consequences on the interacting Chern number of 
the four-dimensional quantum Hall effect. In section [TTT| we study the effect of frequency domain winding number on 
the Z2 index of the three-dimensional topological insulators. We discuss issues on the practical calculations of the 
interacting Green's function in section |IV[ 



II. LOCAL SELF-ENERGY APPROXIMATION AND FREQUENCY DOMAIN WINDING NUMBERS 

We assume the self-energy is local and show how the local fluctuation affects the topological properties. Under 
this assumption the self-energy is independent of momentum k, i.e. E(a;,k) « I](w) but preserves full frequency 
dependencd"^^'. The local self-energy approximation has been widely adopted in the investigations of strongly-correlated 
systems under the name dynamical-mean-field theory (DMFT)^^ and it was proofed to be an accurate approximation 
for high dimensions22l23. 

In the following we show that with nontrivial frequency dependent self-energy, the frequency-domain winding 
number (FDWN) could affect the topological class of the system. Specifically, the interaction induced fluctuation 
effect could change the topological classes without developing LRO. 

Without loss of generality, we assume the Hamiltonian could be expanded on the Dirac gamma matrices iJk = 
Sa=i ^k-'^"' ''■k is a five component vector whose module is denoted as |/ik|. We also introduce the normalized 
vector /i^ = /ik/|/ik|- Under the local self-energy assumption the Green's function could be written as G^^ = 
iw — /i^r° — Tj(liS) = Gatom ~ ^'\^°' i whcrc we introduced the atomic Green's function G~^^^ = ito — S(w) which 
contains all of the frequency dependence. Several lines of algebra^^ shows that the interacting Chern number is 

n=\j d^k r p^ .^i^^^^r .... eabcdehldkXdk^KdkXdk.K (2) 

Denote z = Gatom (^) ^^^^ notice that it defines a mapping from real axis (frequency tu) to the complex plane, we 
could finish the frequency integration 

dw 9<^Gatom r'°°'' dz 



2^» (Gatom - \hu\'r A(-oo) 2m (Z2 _ \h^\2)3 

dz 1 



c 27ri (z2 _ \h^,\2)3 
^ 7 (3) 
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In the second step we appended appropriate large circles to render the integration on a closed contour C, see FigjT] 
Since G~tom is parameterized by w, the contour could be thought of as the path of Gatom ('^) o'^ t^^ complex plane 
as w varies from — cx) to cxd. The result of the contour integration is determined by the singularities that contour C 
encloses. Notice that the residues at the two singularities =F|^k| are zfc -^g,^ ,5 respectively, we introduce 7 to denote 
their winding number differences^^. Plug the result Eqnlslback into Eqnl2l we have 

3 /■ ., » „ ^ ^ 

n = 7X—^ d'^keabcdeh'^dk^hidkh^dk^hldkX (4) 

Eqn. |4] is the central result of the paper. FDWN 7 describes the winding of the atomic Green's function on the 
complex plane. It appears as a pre-factor of the momentum space TKNN number ^'^. In other words, the interacting 
Chern number n is determined by the winding numbers in both the frequency and momentum spaces. Noninteracting 
case has 7=1 since the frequency dependency of the noninteracting Green's function comes trivially from ico, see 
FigjITa). The presence of interacting induced self-energy would arguably give more complex contours and thus change 
the value of FDWN, see Fig|lTb-c). By this mechanism, local fluctuation effect changes the topological class through 
a multiplicity integer factor. 
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FIG. 1: (a-c): The mapping uj t-^ G^ 



appended with appropriate large circles defines a contour in the complex plane. 



The result of the contour integration gives FDWN. Cross denotes the singular points of the integrand 
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means the contour wraps around twice, (d-g): The imaginary part of the PHS (d) G. 
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They have FDWNs 1, 0, 2, 1 respectively. The result could be read off from the number of right and left 



moving intersections of the curve with a horizontal line: 7 equals to the difference of them. One could also get the results from 
the complex plane contour integrations with contours drawn in (a),(b),(c),(a) respectively. 



In case of the particle-hole symmetric (PHS) case, the meaning of FDWN is even simpler and appealing. In fact, 



with PHS the self-energy is purely imaginary and so does the atomic Green's function, i.e. 
real function. We redo the frequency integration in EqnJ2] 
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if{uj) with / a 



dw S^Gg^to^ 
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(5) 



The contribution purely comes from the arctan functions. Notice that u 1— > /(w) defines a mapping from M to M, we 
could view ±(X) as the identical points and compactification the space M into a circle S^ with stereographic projection. 
Then the function / maps S^ to S^ and one has well defined winding numbers deg{/} of the mapping. One nice thing 
at the PHS case is the FDWN could be calculated by simply countering the right moving and left moving intersection 
of /(w) with a given a horizontal line, see Figll] For example, in the noninteracting case, /(cj) = a; indicates 7 = 1, 
T^'-jlla). Nontrivial self-energies like I](w) - ^' ^ ' ^ ^ "■^— --^ ^^"^ - ^ '^ ^ ^ ^^.^u. ^, - o 
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{f ^ (jj — 4f ) results 7 = 2 and 



{f = UJ+ ^fli ) gives 7 = 1. The resulting FDWN are in accodance with the complex integration 



in the complex plane (Figll|a-c)). These results are also been validated by direct numerical integrating 
EqnIT] see fellowing for discussions on the numerical details. 

In general, the interacting Green's functions calculated from numerical solvers may have no analytical expression 
and have much complicated behaviors. The above rules still apply: one can simply read out the FDWN by plotting 
the Green's function and counting the number of intersections or winding numbers. 

We then discuss some physical consequences of the FDWN. First, it indicates that the local fluctuations can affect the 
topological properties of the system by developing nontrivial frequency dependency of the Green's functions. Second, 
through this mechanism, the topological phase could be broken without developing LRO. This phenomenon has been 



uncovered in previous studies on the interacting Haldane modeP. Moreover, T,{uj) = J- is actually corresponds 
to the Mott insulating statea^S. The fact of it has zero FDWN indicates that Mott insulating state in this sense 
has zero Chern number. To clarify potential misconceptions, the topological Mott insulator state refereed ii^ is 
actually a mean-field bond-order wave state with trivial FDWNs 7 = 1. We would also like to remark that in 
many of the proposed interaction induced topological insulating phase d^^"^ ^ interaction effects are mostly considered 
to give nontrivial momentum dependence of the mean-field self-energy Sj., which is effectively absorbed into the 
noninteracting Hamiltonian iJk- 

III. THREE DIMENSIONAL Z2 TOPOLOGICAL INSULATORS 

In this secti on w e show that the FDWN is also relevant to the Z2 index of the three dimensional topological 
insulators. Ref.'^'^'^'show that the three dimensional Z2 topological insulators are closed related to the four dimensional 
Chern insulators through the dimension reduction. Given noninteracting Hamiltonian i/j. (k now denotes three 
dimensional momentum), the three dimensional Z2 index could be calculated through a similar formula as Equfl] We 
extend iJk to four dimension^SI by introducing a pumping parameter kx and H^ — >■ ^k + k\T^ (The choice preserves 
the time reversal symmetry when kx — 00, see^Sl). Under the local self-energy approximation, we integrate over the 
frequency uj and the pumping parameter kx and get. 



n 



3 f f°° 1 



= 7 X £ (6) 

The above formula shows that the FDWN 7 appears as a multiply factor with the momentum space Z2 index £ 
(which is defined upon module 2). As an example, the noninteracting Z2 topological insulators has odd momentum 
space Z2 index multiplied with 7 = 1. 

Notice that to get nontrivial (odd) Z2 index, we need both 7 and £ to be an odd number. This leads to the 
conclusion that the only way FDWN plays a role of changing topological classes is to break down a noninteracting 
Z2 topological insulators, i.e., there is no way to develop a topological nontrivial states from an trivial noninteracting 
state, (since even integer times any integer still gives an even integer numbers) 

It is interesting to look at the three representative cases of the topological trivial states: (a), odd 7 even £, (b). 
even 7 even £ and (c). even 7 odd £. Physically, they correspond to (a), noninteracting band insulators, (b). Mott 
insulators and (c). a state with nontrivial momentum space topology, but overtakes by strong interactions induced 
even FDWN. Notice that in the last case it is not necessary to have LRO. It may be a candidate description of the 
spin liquid state uncovered in the studies of the KMH modeP"Sl. It is interesting to study the FDWN inside the spin 
liquid phase. See the following for discussions on the effect of the spatial fluctuations and practical calculations. 

IV. DISCUSSIONS 

In this section, we comment on the practical calculation of the interacting Green's functions and the techniques for 
the frequency-momentum integration in EqnIT] 

Firstly, under the local assumption, it's suitable to study the interaction effect on topological insulators using 
DMFT. Since the FDWN only depends on the overall shape of the self-energy, simple impurity solvers like IPTZIEII 
will do the trick. And since only the Matsubara Green's function is needed, the advanced CTQMC solvera^SlEil could 
play major role in exact solution of the auxiliary impurity problems. 

Secondly, under the local self-energy approximation, one could get similar results as EqnH] and and Eqnj6] for the 
two dimensional Chern and Z2 insulators. However, one caveat is that the local self-energy assumption may be 
not valid in low dimensions. Nevertheless, the FDWN could still be the driving force for some of the interaction 
induced topological transitions. In case of the strong momentum dependence self-energies, diagrammatic calculations 
such as fluctuation-exchange approximation (FLEX)^^ could be used to calculate the self-energy in the perturbative 
region. Another direction is using the dynamical cluster approximation^ or diagrammatic Monte Carlo methodP^ 



to incorporate the momentum dependence of the self-energies non-perturbatively. These results will tell how much 
the spatial fluctuations affect the results reported here. It's also interesting to see the synergies of these conventional 
many-body numerical methods with the study of topological insulators though the calculation of its interacting Green's 
functions. 

Thirdly, we remark on the numerics for the integration of EqnIT] In cases such as the frequency and momentum 
dependence of the self-energy interwind or I](a;,k) comes from the above mentioned numerical methods, one need to 
perform the five-dimensional integration in EqnjT] numerically. Since usually the contributions of the integrand are 
concentrated around several particular points in the frequency-momentum space (w = and several special k-points) , 
it advised to use Monte Carlo technique to sample the integrand. In our case, we find the VEGAS algorithrrPS (which 
has been used for more than thirty years in the high-energy physics for calculating multi-dimensional quadrature) is 
particularly suitable for this kind of integration. 

Finally, for realistic topological insulating materials which are interacting, it maybe possible to read out the 
frequency-momentum dependency of the self-energy from the ARPES data and plug into Eqn|l] to tell whether 
the the material is in the topological phase or not. 

V. CONCLUSIONS 

In conclusion, we studied the interaction effect on topological insulators through their single-particle Green's func- 
tions. Under the local self-energy approximation, we uncover a mechanism where the winding in the frequency domain 
changes the topological classes. FDWN provide a mechanism of breaking down the topological phases without develop- 
ing LRO. The finding is in accordance with the ongoing numerical studies on the interacting topological insulators'^ S'. 
We urge the interpretation of the featureless intermediate phases uncovered in these studies based on the FDWN. 
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